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Tidal effects arise from differential and inelastic deformation of a planet by a perturbing 
body. The continuous action of tides modify the rotation of the planet together with its orbit 
until an equilibrium situation is reached. It is often believed that synchronous motion is the most 
probable outcome of the tidal evolution process, since synchronous rotation is observed for the 
majority of the satellites in the Solar System. However, in the XIX th century, Schiaparelli also 
assumed synchronous motion for the rotations of Mercury and Venus, and was later shown to 
be wrong. Rather, for planets in eccentric orbits synchronous rotation is very unlikely. The 
rotation period and axial tilt of exoplanets is still unknown, but a large number of planets have 
been detected close to the parent star and should have evolved to a final equilibrium situation. 
Therefore, based on the Solar System well studied cases, we can make some predictions for 
exoplanets. Here we describe in detail the main tidal effects that modify the secular evolution 
of the spin and the orbit of a planet. We then apply our knowledge acquired from Solar 
System situations to exoplanet cases. In particular, we will focus on two classes of planets, 
"Hot-Jupiters" (fluid) and "Super-Earths" (rocky with atmosphere). 



1. INTRODUCTION 

The occurrence, on most open ocean coasts, of high sea 
tide at about the time of Moon's passage across the merid- 
ian, early prompted the idea that Earth's satellite exerts an 
attraction on the water. The occurrence of a second high 
tide when the Moon is on the opposite meridian was a great 
puzzle, but the correct explanation of the tidal phenomena 
was given by Newton in "Philosophic Naturalis Principia 
Mathematical . Tides are a consequence of the lunar and 
solar gravitational forces acting in accordance with laws of 
mechanics. Newton realized that the tidal forces also must 
affect the atmosphere, but he assumed that the atmospheric 
tides would be too small to be detected, because changes 
in weather would introduce large irregular variations upon 
barometric measurements. 

However, the semi-diurnal oscillations of the atmo- 
spheric surface pressure has proven to be one of the most 
regular of all meteorological phenomena. It is readily de- 
tecta ble by harmonic analysis at a ny station over the world 
(e.g. lChapman and Lindzenl ll970). The main difference in 
respect to ocean tides is that atmospheric tides follow the 
Sun and not the Moon, as the atmosphere is essentially ex- 
cited by the Solar heat. Even though tides of gravitational 
origin are present in the atmosphere, the thermal tides are 
more important as the pressure variations on the ground 
are more sensitive to the temperature gradients than to the 
gravitational ones. 

The inner planets of the Solar System as well as the 
majority of the main satellites present today a spin dif- 
ferent from what is believed to have been the initial one 



(e.g. Goldreich and Sote rll96dlGoldreich and Pealell 968). 
Planets and satellites are supposed to rotate much faster in 
the beginning and any orientation of the s pin axis may be 



allow ed (e.g. iDones and Tremaind Il993t iKokubo and Ida 
2007). However, tidal dissipation within the internal lay- 



ers give rise to secular evolution of planetary spins and or- 
bits. In the case of the satellites, spin and orbital evolution 
is mainly driven by tidal interactions with the central planet, 
whereas for the inner planets the main source of tidal dis- 
sipation is the Sun (in the case of the Earth, tides raised by 
the Moon are also important). 

Orbital and spin evolution cannot be dissociated because 
the total angular momentum must be conserved. As a con- 
sequence a reduction in the rotation rate of a body implies 
an increment of the orbit semi-major axis and vice-versa. 
For instance, the Ea rth's rotation p eriod is increasing about 
2 ms/century (e.g. lWilliamslll990l) . and the M oon is cons e- 



quent ly moving away about 3.8 cm/year (e.g. lDickev et al 
1994). On the other hand, Neptune's moon, Triton, and 



the Martian moon, Phobos, are spiraling down into the 
planet, clearly indicating that the present orbits are not pri- 
mordial, and may have undergone a long evolving process 
from a previous capture from an heliocentric orb it (e.g. 
Mignard[l9^lGoldreicheTal1ll989l:lc"orreiall2009h . Both 



the Earth's Moon and Pluto's moon, Charon, have a sig- 
nificant fraction of the mass of their systems, and there- 
fore they could be classified as double-planets rather than 
as satellites. The proto-planetary disk is unlikely to pro- 
duce double-planet systems, whose origin seems to be due 
to a catastrophic impact of the initial planet with a bod 



of comparable dimensions (e.g. lCanup and Asph aug 2001 



■ay 

31; 



l 



Canupl 2005 ). The resulting orbits after collision are most 
likely eccentric, but the present orbits are almost circular 
suggesting that tidal evolution subsequently occurred. 

The ultimate stage for tidal evolution corresponds to 
the synchronous rotation, a configuration where the ro- 
tation rate coincides with the orbital mean motion, since 
synchronous equilibrium corresponds to the minimum of 
dissipation of energy. However, when the eccentricity 
is different from zero some other configurations are pos- 
sible, such as t he 3/2 spin-orbit resonance observed for 
planet Mercury JColomboll965llGoldreich and PealeT l966; 



ICorreia and Laskarl|2004l) or the chaotic rotation of Hype- 
rion ( Wisdom et al.l 1 1 9841) . When a dense atmosphere is 
present, thermal atmospheric tides may also counterbalance 
the gravitational tidal effect and non-resonant equilibrium 
configurations are also possib le, as it is illustrated by t he ret- 
rograde rotation of Venus (ICorreia and Laskarll200ll) . Ad- 
ditional effects may also contribute to the final evolution 
of the spin, such as planetary perturbations or core-mantle 
friction. 

Despite the proximity of Mercury and Venus to the 
Earth, the determination of their rotational periods has only 
been achieved in the second half of the XX th century, 
whe n it became possible to use radar ranging on the plan- 
ets (IPettengill and Dyed IT965b iGoldsteinl 1 1 964t ICarpenterl 



1964). We thus do not expect that it will be easy to observe 



the rotation of the recently discovered exoplanets. Never- 
theless, many of the exoplanets are close to their host star, 
and we can assume that exoplanets' spin and orbit have 
already undergone enough dissipation and evolved into a 
final equilibrium possibility. An identical assumption has 
been done before by Schiaparelli for Mercury and Venus 
1889b . who m ade pr edictions for their rotations based on 



Darwin' work (1880). Schiaparelli's predictions were later 
proved to be wrong, but were nevertheless much closer to 
the true rotation periods than most values derived form ob- 
servations in the two previous centuries. As Schiaparelli, 
we may dare to establish predictions for the rotation peri- 
ods of some already known exoplanets. We hope that the 
additional knowledge that we gained from a better under- 
standing of the rotation of Mercury and Venus will help 
us to be at least as close to the reality as Schiaparelli was. 
Indeed, observations also show that many of the exoplan- 
ets have very eccentric orbits. In some cases eccentricities 



larger than 0.9 are found (e.g. lNaef et al. 2001; Jone s et al 



2006; iTamuz et al. 2008), which opens a wide variety for 
final tidal equilibrium positions, different from what we ob- 
serve around the Sun. 

In this Chapter we will describe the tidal effects that 
modify the secular evolution of the spin and orbit of a 
planet. We then apply our knowledge acquired from So- 
lar System situations to exoplanet cases. In particular, we 
will focus on two classes of planets, "Hot-Jupiters" (fluid) 
and "Super-Earths" (rocky), which are close to the star and 
therefore more susceptible of being arrived in a final equi- 
librium situation. 
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Fig. 1 . — Andoyer's canonical variables. L is the projection of 
the total rotational angular momentum vector L on the principal 
axis of inertia K, and X the projection of the angular momentum 
vector on the normal to the orbit (or ecliptic) k. The angle between 
the equinox of date 7 and a fixed point of the equator A is the hour 
angle 9, and ifi — 7 N + N 70 is the general precession angle. The 
direction of 70 is on a fixed plane E c0 , while 7 is on the mean 
orbital (or ecliptic) E ct of date t. 



2. MODEL DESCRIPTION 

We will first omit the tidal effects, and describe the spin 
motion of the planet in a conservative framework. The mo- 
tion equatio ns will be obtain ed from an Hamiltonian for- 
malism (e.g. lGoldsteinl ll950) of the total gravitational en- 
ergy of the planet (Sect. 2.1). Gravitational tides (Sect. 2.2) 
and thermal atmospheric tides (Sect. 2.3) will be described 
later. We also discuss the impact of spin-orbit resonances 
(Sect. 2.4) and planetary perturbations (Sect. 2.5). 

2.1. Conservative motion 

The planet is considered here as a rigid body with mass 
m and moments of inertia A < B < C, supported by the 
reference frame (I, J, K), fixed with respect to the planet's 
figure. Let L be the total rotational angular momentum and 
(i, j, k) a reference frame linked to the orbital plane (where 
k is the normal to this plane). As we are interested in the 
long-term behavior of the spin axis, we merge the axis of 
figure K with direction of the angular momentum L. In- 
deed, the average of K coincides with L/ L up to J 2 , where 
cos J = L • K dBoue and L askar 2006). J is extremely 
small for large rocky planets ( J « 7 x 10~ 7 for the Earth), 
being even smaller for Jupiter-like planets that behave as 
fluids. The angle between K and k is the obliquity, e, and 
thus, cos e = k • K (Fig.[T|i. 

The Hamiltonian of the motion can be written using 
canonical Andoyer's act ion variables (L , X) and their con - 
jugate angles {9, —tp) dAndoyerl 1 19231: iKinoshital Il977l) . 
L = L ■ K = Cuj is the projection of the angular mo- 



2 



mentum on the G axis, with rotation rate to = 8 — ip cos e, 
and X = L ■ k is the projection of the angular momentum 
on the normal to the ecliptic; 8 is the hour angle between 
the equinox of date and a fixed point of the equator, and 
if) is the general precession angle, an angle that simultane- 
ously accounts for the precession of the spin axis and the 
orbit (Fig.Q). 

2.1.1. Gravitational potential 

The gravitational potential V (energy per unit mass) gen- 
erated by the planet at a generic point of the space r, ex- 
panded in degree tw o of R/r, where R is the pla net's ra- 
dius, is given by (e.g. lTisserandl[ l89lt ISmarlfl95l : 



V(r) 



Gm 
r 



G(B - A) 



G(C - A) 



P 2 (f-J) 



P 2 (f-K) 



(1) 



where r = r/r, G is the gravitational constant, and 
P^(x) = (3a; 2 — 1) /2 are the Legendre polynomials of de- 
gree two. The potential energy U when orbiting a central 
star of mass to* is then: 



U = m*V(r) . 



(2) 



For a planet evolving in a non-perturbed keplerian orbit, 
we write: 



f = cos(ti7 + v)i + sin(tz7 + v)j 



(3) 



where w is the longitude of the periapse and v the true 
anomaly (see Chapter 2: Keplerian Orbits and Dynamics). 
Thus, transforming the body equatorial frame (I, J, K) into 
the orbital frame (i, j, k), we obtain (Fig . [0 : 



r ■ J = — cos w sin ( 



sin w cos a cos e 



f K 



(4) 



sin w sm e , 



where w = m + tp + vis the true longitude of date. The 
expression fo r the potential energy (Eq.|2]i becomes (e.g. 
Correiall2006l) : 



U = — 



Gram* GCm* 

1 3 — EdPi (sin w sin e) 



3Gm* , j 
-^^(B-A)F(8,w,e) 
8r 



(5) 



where 



F{6, w,e) = 2 cos(26» - 2w) cos 4 - 



-2 cos(26» + 2w) sin 4 ( - ) + cos(26») sin 2 e , (6) 



and 



E d - ' ;i 



G 
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where Ed is the dynamical ellipticity, and fc/ is the fluid 
Love number (pertaining to a perfectly fluid body with the 
same mass distribution as the actual planet). The first part of 
Ed (Eq. [7]) corresponds to the flattening in hydrostatic equi- 
librium EambeckH l980). and 5 Ed to the departure from this 
equilibrium. 

2.1.2. Averaged potential 

Since we are only interested in the study of the long- 
term motion, we will average the potential energy U over 
the rotation angle 8 and the mean anomaly M: 



U = 



1 

i^ 2 



U dM d8 



(8) 



However, when the rotation frequency uj s» 8 and the mean 
motion n = M are close to resonance (lj pn, for a 
semi-integer value p), the terms with argument 2(8 — pM) 
vary slowly and must be retained in the expansions (e.g. 
Murr ay and Dermott|[l999l) 



cos(2<9) 



1 



G{p, e) cos 2(8 -pM) , (9) 



and 



cos(26< - 2w) 



+ 00 



= V H(p,e) cos 2(8 - P M) , (10) 
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where a and e are the semi-major axis and the eccentricity 
of the planet's orbit, respectively. The functions G(p, e) and 
H (p, e) can be expressed in power series in e (TableQ}. The 
averaged non-constant part of the potential U becomes: 
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— = — a- 
C 



uj x 
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where x = X/L = cos e, <f> — w + ip, 
3Gto* Ed 3 n 2 



2a 3 (l-e 2 ) 3 / 2 w ~ 2 u 
is the "precession constant" and 

3Gm* B — A 3 



(1 



- 3 ' 2 E d (12) 
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(13) 



For non-resonant motion, that is, when (B — A)/C s» 
(e.g. gaseous planets) or |w| 3> pn, we can simplify expres- 
sion ( fTTT ) as: 



U 
G 



(14) 
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Table 1 

Coefficients of G(p, e) and H (p, e) to e 4 . 
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The exact expression of the coefficients is given by G(p, e) = ^ (^) exp(i2pM) dM and 
H (p, e) = i /* (f ) 3 exp(i 2i/) exp(i 2pM) dM. 



2.7.3. Equations of motion 

The Andoyer variables (L, 6 Q and (X, — ijj) are canon- 
ical^ conjugated and thus (e.g. Goldst ein! 1950t Kinoshital 
1977) 



dL 
~dt 



dU 
~d0 



dX 
~~dt 



dU (hp_ 
~£hp ' ~dt 



dU 
dX 



(15) 



Andoyer's variables do not give a clear view of the spin 
variations, despite their practical use. Since to = L/C and 
cos e = x — X/L the spin variations can be obtained as: 
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(17) 



For non-resonant motion, we get from equation ( fT4l : 

dw dx dip 

— = — = and —=ai. (18) 

dt dt at 

The spin motion reduces to the precession of the spin vector 
about the normal to the orbital plane with rate a x. 



2.2. Gravitational tides 

Gravitational tides arise from differential and inelastic 
deformations of the planet due to the gravitational effect of 
a perturbing body (that can be the central star or a satel- 
lite). Tidal contributions to the planet evolution are based 
on a very general formulation of the tidal potential, initi- 
ated by George H. Darwin (1880). The attraction of a body 
with mass m+ at a distance r from the center of mass of the 
planet can be expressed as the gradient of a scalar pot ential 
V, which is a sum of Legend re polynomials (e.g. iKaula 
1964; Efroimskv and Williams 2009): 



V=5> = -^£(7) fl(cosS), (19) 

1=2 1=2 ^ ' 

where r' is the radial distance from the planet's center, and 
S the angle between r and r'. The distortion of the planet 
by the potential V' gives rise to a tidal potential, 



(20) 



1=2 



where (V g )i — kiV{ at the planet's surface and ki is the 
Love number for potential (Fig.|2]i. Typically, k 2 ~ 0.25 
for Earth-lik e planets, and k% ~ 0.40 for giant planets 
dYoderlll995l) . Since the tidal potential (V g )i is an ^ th de- 
gree harmonic, it is a solution of a Dirichlet problem, and 
exterior to the planet it must be proportional to (e.g. 
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Abramowitz and Stegunlll972t lLambeck 1980). Further- 
more, as upon the surface r' = R <C r, we can retain in 
expression (l20l only the first term, I = 2: 



V a 



Gm± 
' R 



P 2 (cosS). (21) 



In general, imperfect elasticity will cause th e phase angle 
of V g to lag behind that of V dKaulal [l964h by an angle 



S g (a) such that: 



2S g (a) = aAt g (a) , 



(22) 



At g (a) being the time lag associated to the tidal frequency 
(7 (a linear combination of the inertial rotation rate u> and 
the mean orbital motion n) (Fig. [3]). 

2.2.1. Equations of motion 

Expressing the tidal potential given by expression (fJTJ 
in terms of Andoyer angles (6, ijj), we can obtain the con- 
tribution to the spin evolution from expressions ( fT3T > using 
Ug = m'Vg at the place of U, where m! is the mass of the 
interacting body. As we are interested here in the study of 
the secular evolution of the spin, we also average U g over 
the periods of mean anomaly and longitude of the periapse 
of the orbit. When the interacting body is the same as the 
perturbing body (to/ = to*), we obtain: 



ds 
dt 



duo GmlR 5 
~dt ~ Ca 6 

GmlR 5 sine 
~ Cc? uT 
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(23) 



^b g (o-)£*(x,e) , (24) 



where the coefficients Sl%(x, e) and ££(x, e ) stre polynomi- 
als in the eccentricity ( Kaulal[l964 ). When the eccentricity 
is small, we can neglect the terms in e 2 , and we have: 
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+ 
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(25) 
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6 r ( W )|x 3 



2n)§(l 
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b T (2uj - 2n)4(l +x) 3 
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x) 
x) 



(26) 

The coefficients b T (cr) are related to the dissipation of 
the mechanical energy of tides in the planet's interior, re- 
sponsible for the time delay At g (<r) between the position 



of "maximal tide" and the sub-stellar point. They are re- 
lated to the phase lag S g (a) as: 

b g (a) = ki sin 28 g (a) — k2Sm(aAt g (a)) , (27) 

where r = g for gravitational tides. Dissipation equations 
(l23l and (1241) must be invariant under the change (lu, x) by 
(— u>, — x) which imposes that b(a) = —b(—a), that is, b(a) 
is an odd function of u. Although mathematically equiva- 
lent, the couples (w,x) and (—ui, —x) correspond to two 



different physical situations (ICorreia and Las kar 2001). 

The tidal potential given by expression OTl i can also 
be directly used to compute the orbital evolution due to 
tides. Indeed, it can be seen as a perturbation of the grav- 
itational potential (Eq.Q]), and the contributions to the or- 
bit are computed using Lagrange Planetar y equations (e.g. 
Brouwer and Clemencdll96lllKaulall 19641) : 
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(28) 



V dM dm 



(29) 



where /i = mm^j (to + to*) to is the reduced mass. 
We then find for the orbital evolution of the planet: 
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(31) 



where the coefficients A%{x,e) and E§.(x,e) are again 
polynomials in the eccentricity. When the eccentricity is 
small, we can neglect the terms in e 2 , and we have: 



& T (2n)^(l-x 2 ) 2 
6 T (w-2n)|(l-x 2 )(l + x) 2 
6 T (w + 2n)|(l -x 2 )(l -x) 2 
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(32) 
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Fig. 2. — Gravitational tides. The difference between the gravitational force exerted by the mass m on a point of the surface and the 
center of mass is schematized by the arrows. The planet will deform following the equipotential of all present forces. 




2.2.2. Dissipation models 

The dissipation of the mechanical energy of tides in the 
planet's interior is responsible for the phase lags 8(a). A 
commonly used dim ensionless measure of tidal d amping is 
the quality factor Q (iMunk and MacDo nald 1960), defined 
as the inverse of the "specific" dissipation and related to the 
phase lags by 



Q(?) 



2irE 

~Ke 



COt2<5(cr) , 



(34) 



where E is the total tidal energy stored in the planet, and 
AE the energy dissipated per cycle. We can rewrite expres- 
sion d2~7T i as: 



k 2 sign(cr) 



sign(cr) 



(35) 



The present Q value for the planets in the Solar system can 
be estimated from orbital measurements, but as rheology of 
the planets is badly known, the exact dependence of b T (a) 
on the tidal frequency a is unknown. Many different au- 
thors have studied the problem and several models for b T (a) 
have been developed so far, from the simplest ones to the 
more complex (for a review see Efroimsk v and Will iams 
2009). The huge problem in validating one model better 
than the others is the difficulty to compare the theoretical 
results with the observations, as the effect of tides are very 



small and can only be detected efficiently after long peri- 
ods of time. Therefore, here we will only describe a few 
simplified models that are commonly used: 



The visco-elastic model Darwinl ( 1908 ) assumed that the 
planet behaves like a Maxwell solid, that is, the planet re- 
sponds to stresses like a massless, damped harmonic oscil- 
lator. It is characterized by a rigidity (or shear modulus) 
p e and by a viscosity v e . A Maxwell solid behaves like 
an elastic solid over short time scales, but flows like a fluid 
over long periods of time. This behavior is also known as 
elasticoviscosity. For a constant density p, we have: 



(36) 



where kf is the fluid Love number (Eq.|7]). r Q = v e /p e 
and Tb — T a (l + WfieR/SGmp) are time constants for the 
damping of gravitational tides. 

The visco-elastic model is a realistic approximation 
of the planet's defo rmation with the tidal frequency (e.g. 
Escribano et al.ll2008l) . However, when replacing expres- 
sion (l36*l l into the dynamical equations (|23T > and d24l i we 
get an infinite sum of terms, which is not practical. As a 
consequence, simplified versions of the visco-elastic model 
for specific values of the tidal frequency a are often used. 
For instance, when a is small, (r^ a) 2 can be neglected in 
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expression (f36b and b g (a) becomes proportional to a. 

The viscous or linear model In the viscous model, it is 
assumed that the response time delay to the perturbation is 
independent of the tidal frequency, i.e., the position of the 
"maximal tide" is shift ed from the sub-stell ar point by a 
constant time lag At g dMignardlll979[ Il980l) . As usually 



we have aAt q <C 1, the viscous model becomes linear: 



fc 2 sin(erAi g ) w k 2 <?At g 



(37) 



The viscous model is a particular case of the visco-elastic 
model and is specially adapted to describe the behavior of 
planets in slow rotating regimes (w ~ n). 

The constant-Q model Since for the Earth, Q changes by 
less than an order of magnitude between the Chandler wob- 
ble period (ab out 440 days) and seism ic periods of a few 
seconds (iMunk and MacDonaldl ll960). it is also common 
to treat the specific dissipation as independent of frequency. 
Thus, 



b g {(i) « sign(cr)fc 2 /Q 



(38) 



The constant- Q model can be used for periods of time 
where the tidal frequency does not change much, as is the 
case for fast rotating planets. However, for long-term evo- 
lutions and slow rotating planets, the constant-Q model is 
not appropriate as it gives rise to discontinuities for a = 0. 

2.2.3. Consequences to the spin 

Although both linear and constant models have some 
limitations, for simplicity reasons they are the most widely 
used in literature. The linear model has nevertheless an im- 
portant advantage over the constant model: it is appropriate 
do describe the behavior of the planet near the equilibrium 
positions, since the linear model closely follows the realis- 
tic visco-elastic model for slow rotation rates. The equa- 
tions of motion can also be expressed in an elegant way, 
so we will adopt the viscous model for the remaining of 
this Chapter, without loss of generality concerning the main 
consequences of tidal effects. 

Using the approximation d37b in e xpressions d23l and 
d24l) . we simplify the spin equations as (iCorreia and Laskat 
2010, Appendix B): 



and 



where 



K n I 1 
~C /l(e) " 2 



cos e iij \ 

--/ 2 e)cose , (39) 
n / 



Kn 
e ~ — — sin e 

Cuj 



cose- -Me) 



A(e) 



1 + 3e 2 + 3e 4 / 



h(e) = 



(l_ e 2)9/2 ' 

l + 15e 2 /2 + 45e 4 /8 + 5e 6 /16 



2nG 



(l-e a ) 



(40) 
(41) 
(42) 
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Fig. 4. — Evolution of the equilibrium rotation rate u) e /n — 
/2(e) //i(e) with the eccentricity when e = 0° using the viscous 
model (Eq.l44t. As the eccentricity increases, ui £ also increases. 
The gravitational tides lead the planet to exact resonance when the 
eccentricity is respectively e 1 / 1 — 0, e 3 / 2 — 0.284926803 and 
e 2 /i = 0.392363112. 



and 



K = At 



3k 2 GmlR 5 



(43) 



Because of the factor 1/w in the magnitude of the obliq- 
uity variations (Eq.l40li. for an initial fast rotating planet the 
time-scale for the obliquity evolution will be longer than the 
time-scale for the rotation rate evolution (Eq.[39l. As a con- 
sequence, it is to be expected that the rotation rate reaches 
an equilibrium value earlier than the obliquity. For a given 
obliquity and eccentricity, the equilibrium rotation rate, ob- 
tained when ill = 0, is then attained for (Fig. [4]): 



n 



/2(e) 2 cose 



fi (e) 1 + cos 2 £ 



(44) 



Replacing the previous equation in the expression for obliq- 
uity variations (Eq.l40"l>. we find: 



Kn 



sine 



(45) 



We then conclude that the obliquity can only decrease by 
tidal effect, since i < 0, and the final obliquity always tends 
to zero. 

2.2.4. Consequences to the orbit 

As for the spin, the semi-major axis and the eccentricity 
evolution can be obtained using the ap proximation (f37T > in 
expressions d32b and d33T >. respectively jCorre ia 2009): 



— (/a (e) cose- - /3(e)) 



and 

where 

/s(e) = 



fin 
9K /ll 



fj,a 2 \18 
l + 31e 2 /2 + 255e 4 



/ 4 (e)cose- - / 5 (e) ) e, 
n 



(46) 



(47) 



185e 6 /16 + 25e 8 /64 



(l_ e 2)15/2 



(48) 
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hie) 



. . . 1 + 3e 2 /2 + e 4 /8 

/4(6) = (l- e 2)5 < 

1 + 15e 2 /4 + 15e 4 /8 + 5e 6 /64 
(1-e 2 ) 13 / 2 



(49) 



(50) 



The ratio between orbital and spin evolution time-scales 
is roughly given by C /(pa 2 ) <C 1, meaning that the spin 
achieves an equilibrium position much faster than the orbit. 
Replacing the equilibrium rotation rate (Eq.l44b with e = 
(for simplicity) in equations d46l ) and ( [47) , gives: 



7K 
pa 



/e(e)e 2 



IK 

2pa 2 



/ 6 (e)(l-e 2 )e, 



(51) 



(52) 



where / 6 (e) = (1 + 45e 2 /14 + 8e 4 + 685e 6 /224 + 
255e 8 /448 + 25e 10 /1792)(l~e 2 )- 15 / 2 /(l+3e 2 + 3e 4 /8) . 
Thus, we always have a < and e < 0, and the final ec- 
centricity is zero. Another consequence is that the quan- 
tity a(l — e 2 ) is conserved (Eq. llOU . The final equilibrium 
semi-major axis is then given by 



a/ = a(l — e 2 ) , 



(53) 



which is a natural consequence of the orbital angular mo- 
mentum conservation (since the rotational angular momen- 
tum of the planet is much smaller). Notice, however, that 
once the equilibrium semi-major axis a/ is attained, the 
tidal effects on the star cannot be neglected, and they govern 
the future evolution of the planet's orbit. 

2.3. Thermal atmospheric tides 

The differential absorption of the Solar heat by the 
planet's atmosphere gives rise to local variations of tem- 
perature and consequently to pressure gradients. The mass 
of the atmosphere is then permanently redistributed, ad- 
justing for an equilibrium position. More precisely, the 
particles of the atmosphere move from the high tempera- 
ture zone (at the sub-stellar point) to the low temperature 
areas. Indeed, observations on Earth show that the pressure 
redistributi on is essentially a superposi tion of two pressure 
waves (see IChapman and Lindzen 1970): a daily (or diur- 
nal) tide of small amplitude (the pressure is minimal at the 
sub-stellar point and maximal at the antipode) and a strong 
half-daily (semi-diurnal) tide (the pressure is minimal at the 
sub-stellar point and at the antipode) (Fig.|5J. 

The gravitational potential generated by all of the parti- 
cles in the atmosphere at a generic point of the space r is 
given by: 

V. = -G[*", (54) 
J \r-r'\ 

(M) 

where r' = (r',6',cp') is the position of the atmosphere 
mass element dM with density p n (r') and 



Assuming that the radius of the planet is constant and 
that the height of the atmosphere can be neglected, we ap- 
proximate expression d55l l as: 

i? 2 

dM = — p s (6', ip', t) smO' dff dtp' , (56) 
9 

where g is the mean surface gravity acceleration, and p s the 
surface pressure, which depends on the stellar insolation. 
Thus, p s depends on S, the angle between the direction of 
the Sun and the normal to the surface: 

+00 

Ps {6',ip')= Ps {S)=Y,PlPi(™sS) , (57) 
1=0 

where Pi are the Legendre polynomials of order I and p\ 
its coefficients. Developing also \r — r'| _1 in Legendre 
polynomials we rewrite expression d54b as: 



V a 



lg? 3 
21 + 1 



Pi 



R 



1+1 



Pi(cosS), (58) 



where p is the mean density of the planet. Since we are 
only interested in pressure oscillations, we must subtract 
the term of constant pressure (I = 0) in order to obtain the 
tidal potential. We also eliminate the diurnal terms (I = 
1) because they correspond to a displacement of the center 
of mass of the atmosphere bulge which has no dynamical 
implications. Thus, since we usually have r ^> R, retaining 
only the semi-diurnal terms (I — 2), we write: 



V a = - 3 - P 4(-) P 2 (cosS) 
bp \ r ' 



(59) 



2.3.1. Equations of motion 

Using the same methodology of previous sections, the 
contributions of thermal atmospheric tides to the spin evo- 
lution are obtained from expressions (TT3T > using U a = m*V a 
at the place of U: 



dco 
~dt 



3m*i? d 



hCpa 



3 /L^i 



b a (a)n a a (x,e) 



de 
dt 



3m*i? sine y-y 



U! 



b a (*)£Z(x,e) 



(60) 



(61) 



where the terms 0°(a;, e) and £%(x, e) are also polynomi- 
als in the eccentricity, but different from their analogs for 
gravitational tides (Eqs.l23landl24b. Nevertheless, when ne- 
glecting the terms in e 2 , they become equal and are given 
by expressions d25l and d26l >. respectively (with r = a). 

For thermal atmospheric tides there is also a delay before 
the response of the atmosphere to the excitation (Fig. [6]). We 
name the time delay At a (a) and the corresponding phase 
angle 6 a (cr) (Eq.l22l). The dissipation factor b a {a) is here 
given by: 



dM = p a (r')r' 2 sin 9' dr' dO' dtp' 



(55) 



b a (a) = p 2 ((r)sm2S a (a) = p 2 (a) sm(aAt a (a)) . (62) 
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Siebert ( 1961 ) and Chapman and Lindzen ( 1970b have 
shown that when 

|p 2 (<7)| <|>o, (63) 

the amplitudes of the pressure variations on the ground are 
given by: 



- / n • 7 - 7 - 1 J<? 
p 2 (cr) = i-po V ■ v a — 



(64) 



where 7 = 7/5 for a perfect gas, v is the velocity of tidal 
winds, J rj the amount of heat absorbed or emitted by a unit 
mass of air per unit time, and Hq is the scale height at the 
surface. We can rewrite expression (l64l as, 



7 - 
r-rPo 
M 



V-t>„ 



7-1 J* 
7 .9#o 



e ±if 



|p 2 (cr)|e 



where the factor e ±; a can be seen as a supplementary phase 
lag of ±tt/2: 

&«,(*) = |P2(cr)|sin2^ a (a)±|) 

= -\p 2 (a)\sm25 a (a) . (66) 

The minus sign above causes pressure variations to lead the 
Sun whenever 6 a (a) < tt/ 2 ( Chapman and Lindzen 1970t 



Dobrovolskis and Ingersolll 1980h rFig.l6b. 



2.3.2. Dissipation models 

Unfortunately, our knowledge of the atmosphere re- 
sponse to thermal excitation is still very incomplete. As 
for the gravitational tid es, models are developed t o deal 
with the unknown part. iDobrovolskis and Ingersoll (1980) 
adopted a model called "heating at the ground", where they 
suppose that all the stellar flux absorbed by the ground F s is 
immediately deposited in a thin layer of atmosphere at the 
surface. The heating distributing may then be written as a 
delta-function just above the ground: 



J(r) = JLF s 5{r~0+) 
Pa 



(67) 



Neglecting v over the thin heated layer, expression d64l > 
simplifies: 



\M<r) 



5 7 - 1 F s 5 7 gF s 



16 



16 \a\ c p T s 



(68) 



where the factor 5/16 represents the second- degree har- 



moni c component of the insolation contribution (IDobrovolskis and Ingersol 
19801) and c p is the specific heat at constant pressure. 

Nevertheless, according to expression d68l ). if a = 0, 
the amplitude of the pressure variations p 2 {o) becomes in- 
finite. The amplitude cannot grow infinitely, as for a tidal 
frequency equal to zero, a steady distribution is attained, 
and thus p2(0) = 0. Indeed, expression d64i > is not valid 
when a s» because the condition d63l > is no longer ver- 
ified. Using the typically accepted values for the Venu- 
sian atmosphere, c„ » 1000 K kg~ 1 s~ 1 , T s 730 K and 



(65) F s w 100 Wm~ 2 dAvduevskii et al.ll 1 976h we compute: 



\p 2 (a)\ w 10- 4 j5 01 



(69) 



which me ans that for a ~ n, the "heating a t the ground" 
model of Dobrovolskis and Ingersolll ( 1980h . can still be 
applied. Since we are only interested in long-term be- 
haviors we can set P2(er) = whenever \a\ <C n/100. 
Moreover, for tidal frequencies a ~ the dissipation lag 
sin(crAt tt (cr)) w a/S.t a (a) also goes to zero. We expect 
that further s tudies about synchronous exoplanets ' atmo - 
spheres (e.g. Ijoshi et al.l fl997 : Arras and Socrates 2010), 
may provide a more accurate solution for the case a ~ 0. 

In presence of a dense atmosphere, another kind of tides 
can arise: the atmosphere pressure upon the surface gives 
rise to a deformation, a pressure bulge, that will also be 
affected by the stellar torque. At the same time, the atmo- 
sphere itself exerts a torque over the planet's bulges (gravi- 
tational and pressure bulge). Nevertheless, we do not need 
to take into account additional tidal effects as their con- 
sequences upon the dynamical equations can be neglected 
( Hinderer et al.ll 1 987 : Correia and Laskarll2003al) . 



2.4. Spin-orbit resonances 

A spin-orbit resonance occurs when there is a commen- 
surability between the rotation rate uj and the mean motion 
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of the orbit n (Eqs.|9land[T0l). The synchronous rotation of 
the Moon is the most common example. After t he discov- 
ery o f the 3/2 spin-orbit resonance of Mercury dColombo 
1 19651). spin-orbit resonances were studied in great detail 
dColombo and Shapiro! 1 19661 [Goldreich and Peald 1966 ; 



Counselma n and Shapiro! 1 1970t ICorreia and Laskai 



2004 



2009ll2010h ." When resonant motion is present we cannot 
neglect the terms in (3 in expression fTTV Assuming for 
simplicity a low obliquity (x w 1) we obtain a non-zero 
contribution for the rotation rate (Eq.[T6l>: 



dio 
~dt 



-(3 H(p, e) sin 27 , 



(70) 



where 7 = 8—pM—cf). The rotation of the planet will there- 
fore present oscillations around a mean value. The width of 
the corresponding resonance, centered at u> = pn, is: 



Aw = ^2/3H(p,e) 



(71) 



Due to the tidal torque (Eq.l23l. here denoted by T, 
the mean rotation rate does not remain constant and may 
therefore cross and be captured in a spin-orbit resonance. 
Goldreich and Peald (11966b computed a simple estimation 
of the capture probability P cap , and subsequent more de- 
tailed studies proved their expressio n to be essentially cor- 
rect (for a review, see lHenrardl 1993b . Since the tidal torques 



ca n usually be described by me ans of the torques considered 
bv lGoldreich and Peald (1 1966b . we will adopt here the same 
notations. Let 



where Aw is the resonance width (Eq.[7Tl). In the slow rota- 
tion regime (a; ~ n), where the spin encounters spin-orbit 
resonances and capture may occur, we compute for the vis- 
cous tidal model (Eq.[39l: 



P — 2 

J cap — ^ 



P 



2x /2(e) \ nir 



1 + x 2 /i(e) / 2 Aw 



(74) 



2.5. Planetary perturbations 

Like in the Solar System, many exoplanets are not alone 
in their orbits, but belong to multi-planet systems. Because 
of mutual planetary perturbations the orbital parameters of 
the planets do not remain constant and undergo secular vari- 
ations in time (Chapter 10: Non-Keplerian Dynamics). An 
important consequence for the spin of the planets is that 
the reference orbital plane (to which the obliquity and the 
precession were defined) will also present variations. We 
can track the orbital plane variations by the inclination to 
an inertial reference plane, /, and by the longitude of the 
line of nodes, fi. Under the assumption of principal axis ro- 
tation, the ener gy perturbation attached to an inertial frame 
can b e written (Kinoshita ll977 ; Neron de Surgv and Laskai 
1997): 



[X(l — cos I) — L sin e sin / cos <p] 



L sine sin ip — , 
dt 



dn 

dt 



(75) 



T 



-K V 



(72) 



where K and V are positive constant torques, and 7 = 0; — 
pn. The probability of capture into resonance is then given 
by dGoldreich and Pealel ll966b : 



Pc 



1 + nVn/Aoj 1 



(73) 



where ip = — fl — ip. 

Although the Solar System motion is chaotic dLaskar 
19891 1 1990b . the motion can be approximated over several 
million of years by quasi-periodic series. In particular for 
the orbital elements that are i nvolved in the precession driv- 
ing terms (Eq.l75ll. we have ( Laskar and Robutellll993 ): 



dl 
dt 



dn 

dt 



Mil / I e in =J2 Jke 



i(i/fct+0fc) 



(76) 
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and 



(1 - cos I) — = ^ £k cos(v k t + ipk) , (77) 



where v k are secular frequencies of the orbital motion with 
amplitude J k and phase fa, and i = We may then 

rewrite expression ( TTSb as 



pp — 



L ^2 [A cos(v k t + ipk) x 



-JkVl- x 2 sin(v k t + tp + fa) 



■(78) 

Assuming non-resonant motion, from equations ( [TBI and 
dTTI i we get for the spin motion: 



de 
dt 



= ^Jk cos(v k t + ip + fa) 



and 



dip 
dt 



= a cos e — L k cos(u k t - 



cote^ J k sm(v k t + ip + fa) . (80) 

k 



For planetary systems like the Solar System, the mu- 
tual inclinations remain small dLaskari [l990: Correi a et al.l 
2010), and it follows from expressions (|76| > and dTTb that 
the amplitudes of J k and C k are bounded respectively by 



Jk ~ v k I n 



£-k ~ f k I max /2 



(81) 



The term in C k in expression ( f80b for the precession varia- 
tions can then be neglected for small inclinations. 

From expression d79l ) it is clear that a resonance can oc- 
cur whenever the precession frequency is equal to the oppo- 
site of a secular frequency v k (that is, ip — —v k ). Retaining 
only the terms in k, the problem becomes integrable. We 
can search for the equilibrium positions by setting the obliq- 
uity variations equal to zero (e = 0). It follows then from 
expression ( 179b that ip + v k t + fa = ±7r/2, and replacing 
it in expression (f80b with ip = —v k , we get 



(82) 



a cos e sin e + v k sin e ss j7/c cos e 



which gives the equilibrium positions for the spin of 
the planet, generally know n as "Cassini states" (e.g. 
Henrard and Murigandell987l) . Since J k jv k < 1 (Eq.l8TT). 
the equilibrium positions for the obliquity are then: 



tane : 



Jk 



v k ±a 



cose 



Vk 

a 



(83) 



When \a/v k \ <C \a./v\ C rit, the first expression gives states 
2 and 3, while the second expression has no real roots (states 
1 and 4 do not exist). When \a/v k \ ^> \a,/v\ cr u, the first 



expression approximates Cassini states 1 and 3, while the 
second one gives states 2 and 4. States 1, 2 and 3 are sta- 
ble, while state 4 is unstable. Although gravitational tides 
always decrease the obliquity (Eq.l45l>. the ultimate stage 
of the obliquity evolution is to be captured into a Cassini 
resonant state, similarly to the capture of the rotation in a 
spin-orbit resonance (Eq.l73l. 

The complete system (Eqs.l79land[80l) is usually not in- 
tegrable as there are several terms in expression dTol i, but 
we can look individually to the location of each resonance. 
When the resonances are far apart, the motion will behave 
locally as in the integrable case, with just the addition of 
supplementary small oscillations. However, if several res- 
onances overla p, the motion is no longer r egular and be- 
comes chaotic ( Chiri kovll 1 979t lLaskarll 1996b . For instance, 
the present obliquity variations on Mar s are chaotic and can 
vary from zero to nearly sixty degrees dLaskar and Robutel 



(79) ll993l:lTouma and Wisdomfl 993: LaskaretalJ[2004b) 



3. APPLICATION TO THE PLANETS 

The orbital parameters of exoplanets are reasonably well 
determined from radial velocity, transit or astrometry tech- 
niques, but exoplanets' spins remain a mystery. 

The same applies to the primordial spins of the terrestrial 
planets in the Solar System, since very little constraints can 
be derived from the present planetary formation models. In- 
deed, a small number of large impacts at the end of the for- 
mation process will not average, and can change the spin 
direction. The angular velocitie s are also unpredictable, but 



they are usually high, uj ^> n (Dones and Tremaindll993 



AgnoretalJl999l: iKokubo and Idal2007l) . although impacts 
can also form a slow rotating planet (u> ~ n) if the size of 
the typical accreting bodies is much sm aller than the proto- 
planet (e.g. Schlichting and Sarill2007l). The critical angula r 
velocity for rotational instability is lKokubo and I da (2007): 



3.3 



3 gem 



1/2 



hr" 



(84) 



which sets a maximum initial rotation periods of about 
1.4 h, for the inner planets of the Solar System. 

For the Jovian planets in the Solar System no impor- 
tant mechanism susceptible of altering the rotation rate is 
known, but the orientation of the axis may also change by 
secular resonance with the p lanets (e.g. Correia and Laskar 



secular resonance with the planets (e.g. Correia and Laskar 
2003ri IWard and Hamilton! |2004 iBoue and Laskarll2010l) . 
The fact that all Jovian planets rotate fast (Table[3]l seems to 



be in agreement with theoretical predictions (e.g. lTakata and Stevenson 
19961) . 



An empirical relation derived by iMacDonaldl (1964) 
based on the present rotation rates of planets from Mars to 
Neptune (assumed almost unchanged) gives for the initial 
rotation rates 

cj oc m 4/5 i?r 2 . (85) 

Extrapolating for the remaining inner planets, we get initial 
rotation periods of about 18. 9 h, 13. 5 h, and 1 2 . 7 h for Mer- 
cury, Venus and the Earth, respectively, much faster than 
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today's values, which in agreement with the present forma- 
tion theories. 

The above considerations and expressions can also be 
extended to exoplanets. However, since many of the exo- 
planets are close to their host stars, it is believed that the 
spins have undergone significant tidal dissipation and even- 
tually reached some equilibrium positions, as it happens for 
Mercury and Venus in the Solar System. Therefore, in this 
section we will first review the rotation of the terrestrial 
planets, and then look at the already known exoplanets. In 
particular, we will focus our attention in two classes of ex- 
oplanets, the "Hot-Jupiters" (fluid) and the "Super-Earths" 
(rocky with atmosphere), for which tidal effects may play 
an important role in orbital and spin evolution. 

3.1. Solar System examples 

3.1.1. Mercury 

The present spin of Mercury is very peculiar: the planet 
rotates three times around its axis i n the same time as 
it com pletes two orbital revolutions dPettengill and Dyed 
1965). Within a year of the discovery, the stability of 
this 3/2 spin-orbit resonance became understood as the 
result of the solar torque on Mercury's quadrupolar mo- 
ment of inertia combined wi t h an eccentric orbit ( Eq . l70t 
(Col ombo and Shapiro! fl%5 Goldreich and Peale 1966). 
The way the planet evolved into the 3/2 configuration re- 
mained a mystery for long time, but can be explained as the 
result of tidal evolution combined with t he eccentricity vari- 
ations due to planetary perturbations dCorreia and Laskai 



Mercury has no atmosphere, and the spin evolution of 
the planet is therefore controlled by gravitational tidal in- 
teractions with the Sun. Tidal effects drive the final obliq- 
uity of Mercury close to zero (Eq.l45ll. and the averaged 
equation for the rotation motion near the p resonance can 
be written putting expressions d70t and ( 1391 together: 



/2(e) 
/i(e) 



(86) 



where 7 = 9 - pM - <j>, (3' = (3H(p, e), and K' = 
Knfi (e) jC. Note that near the p resonance a contribution 
from core-mantle friction may also be present, but we will 
negle c t it her e (for a full description see lCorreia and Laskai 
2009[ l2Q10h . The tidal equilibrium is achieved when 
duj/dt — 0, that is, for a constant eccentricity e, when 
Lu/n — /2(e)//i(e). In a circular orbit (e = 0) the 
tidal equilibrium coincides with synchronization, while 
the equilibrium rotation rate ui/n = 3/2 is achieved for 
e 3/2 = 0.284927 (Fig.!). 

For the present value of Mercury's eccentricity (e ~ 
0.206), the capture probability in the 3/2 spin-orbit reso- 
nance is only about 7% (Eq.l74l. However, as the eccentric- 
ity of Mercury suffers strong chaotic variations in time due 
to planetary secular perturbations, the eccentricity can vary 
from nearly zero to more than 0.45, and thus reach values 
higher than the critical value e 3 / 2 = 0.284927. Additional 
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Fig. 7. — Simultaneous evolution of the rotation rate (a) and the 
eccentricity (b) of Mercury. In this example, there is no capture 
at the first encounter with the 3/2 resonance (at t w —3.9974 
Gyr). About 100 Myr later, as the mean eccentricity increases, 
additional crossing of the 3/2 resonance occurs, leading to capture 
with damping of the libration dCorreia and Laskaill2006l) . 



capture into resona nce can then occur, at a ny time during 
the planet's history dCorreia and Laskarll2004l) . 

In order to check the past evolution of Mercury's spin, it 
is not possible to use a single orbital solution, as due to the 
chaotic behavior the motion cannot be predicted precisely 
beyond a few tens of millions of years. A statistical study of 
the past evolutions of Mercury's orbit was then performed, 
with the integration of 1000 orbits over 4 Gyr in the past, 
starting with very close initial conditi ons, within the un- 
certainty o f the present determinations dCorreia and Laskai 
I2004ll2009h . 

For each of the 1000 orbital motion of Mercury, the ro- 
tational motion (Eq.l86ll was integrated numerically with 
planetary perturbations. As e(t) is not constant, ui(t) will 
tend towards a limit value Q(t) that is similar to an aver- 
aged value of (/2//i)(e(t)) and capture into resonance can 
now occur more often (Fig. [7]). Globally, only 38.8% of the 
solutions did not end into re sonance, and the final ca pture 
probability distribution was dCorreia and Laskaill2004l) : 



P x/l = 2.2%, P 3/2 = 55.4%, P 2/1 = 3.6% . 

With the consideration of the chaotic evolution of the 
eccentricity of Mercury, the present 3/2 resonant state be- 
comes the most probable outcome for the spin evolution. 
The largest unknown remains the dissipation factor fc 2 Ai 
in the expression of K (Eq.l43l. A stronger dissipation in- 
creases the probability of capture into the 3/2 resonance, as 
uj/n would follow more closely /2(e) / fi(e) (Fig.|7]i, while 
lower dissipation slightly decreases the capture probabil- 
ity. The inclusion of core-mantle friction also increases the 



chances of capture for all resonances (Co rreia and Laskai 
20091 l2010h . 
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Table 2 

Possible final spin states of Venus. 



state 


s 




P (days) 


Ps (days) 


F a + 


0° 


n + lo s 


76.83 


116.75 


F ~ 


0° 


n — u s 


-243.02 


-116.75 


F+ 


180° 


—n — lo s 


-76.83 


116.75 


F- 


180° 


—n + u) s 


243.02 


-116.75 



There are two retrograde states (F and F w ) and t wo prograde states (F^ an d F+). In all cases the synodic period P s is the same 

dCorreia and Laskaj|200lh . 



3.1.2. Venus 

Venus an unique case in the Solar System: it presents a 
slow retrograde rotation, with an obliquity close to 180 de- 
grees and a 243 -day period ( Smithlll963t lGoldsteinlll964 ; 
Carpenter 1970). According to planetary formation theo- 
ries it is highly improbable that the present spin of Venus is 
primordial, since we would expect a lower obliquity and a 
fast rotating planet (Eq.|83Tl. 

The present rotation of Venus is believed to represent 
a steady state resulting from a balance between gravita- 
tional tides, which drives the planet toward synchronous ro- 
tation, and thermally dr iven atmospheric tides , which drives 
the rotation away (e.g. iGold and Soterlll969l) . The conju- 
gated effect of tides and core-mantle friction can tilt Venus 
down during the planet past evolut ion, but requires high val- 
ues o f the initial obliquity (e.g. lDobrovolskisI 19801 Yoder 
19971) . However, the crossing in the past of a large chaotic 
zone for the spin, resulting from s ecular planetary pertur- 
bations dLaskar and Robutellll993l) . can lead Venus to the 
present retrograde configuration f or most initial conditions 
dCorreia and Laskarll200l[|2003bh . 

Venus has a dense atmosphere and the planet is also 
enough close to the Sun to undergo significant tidal dissi- 
pation. Venus' spin evolution is then controlled by tidal 
effects (both gravitational and thermal). Tidal effects com- 
bine d can drive the obliquity either to e = 0° or e = 180° 
dCorreia et al.l 120031) . For the two final obliquity possibili- 
ties, the tidal components become very simplified, with (at 
second order in the planetary eccentricity) a single term of 
tidal frequency a — 2us — 2n for e — and a = 2oj + 2n 
for e — tt (Eq.l25t. Combining expressions d23l and d60l ) 
we can write for the rotation rate: 



sions d23l and d60i >. respectively. Let f(a) be defined as 



duj 
lit 

du> 
~dt 



[K g b g (2uj - 2n) + K a b a (2u - 2n)\ 



[K g b g {2uj + 2n) + K a b a (2cu + 2n)\ 



(87) 

where K g and K a are given by the constant part of expres- 



bg{2a) 
b g (2a) 



(88) 



As b T (a) is an odd function of a (Eq.[25l. f(a) is an even 
function of a of the form Thus, at equilibrium, with 

du) jdt = 0, we obtain an equilibrium condition 



Mu-xn\) = -!$!■ 



(89) 



where x = +1 for e = and x = — 1 for e = tt. Moreover, 
for all commonly used dissipation models / is monotonic 
and decreasing for slow rotation rates. There are thus only 
four possible values for the final rotation rate w/ of Venus, 
given by 



(90) 



Assuming that the present rotation of Venus corresponds 
to a stable retrograde rotation, since u s > the only possi- 
bilities for the present rotation are e = and uj b s = n — u s , 
or e = tt and Lj b s = oj s — n. In both cases, u) a = n + \u> bs\ 
(uj s is thus the synodic frequency). With 



uj obs = 2tt/243.0185 day ; n= 2tt/224.701 day 
we have 

w s = 27r/116.751 day . 



(91) 
(92) 



We can then determine all four final states for Venus (Ta- 
ble O. There are two retrograde states (F _ and F~) and 
two prograde states (Fq~ and F+). The two retrograde 
states correspond to the observed present retrograde state of 
Venus with period 243.02 days, while the two other states 
have a prograde rotation period of 76.83 days. Looking to 
the present rotation state of the planet, it is impossible to 
distinguish between the two states with the same angular 
momentum (Fig. [8). 

In order to obtain a global view of the possible final 
evolutions of Venus' spin, numerical integrations of the 
equations of motion for the dissipative effects (Eqs. l23ll24l 
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Fig. 8. — Final states for a planet with strong atmospheric ther- 
mal tides. The original equilibrium point obtained at synchro- 
nization (uj/n = 1) when considering uniquely the gravitational 
tides, becomes unstable, and bifurcates at e = into two new 
stable fixed points F n ~ and , and at e = tv into F~ and 
dCorreia and Laskaj|200l[|2003bh . 



l60l and IrjTt . with the addition o f planetary perturbations 
(Eqs.| 79landl80ll were performed ( Correia and Laskarl 200 1 . 
l2003bl) . In Figure|9]we show the possible final evolutions 
for a planet starting with an initial period ranging from 3 to 
12 days, with an increment of 0.25 day, and initial obliquity 
from 0° to 180°, with an increment of 2.5 degrees (rotation 
periods faster than 3 days are excluded as they do not allow 
the planet to reach a final rotation state within the age of the 
Solar System). Each color represents one of the possible fi- 
nal states. For high initial obliquities, the spin of Venus al- 
ways evolves into the retrograde final state F~ . It is essen- 
tially the same evolution as without planetary perturbations, 
since none of t he trajectories encounter s a chaotic zone for 
the obliquity dLaskar and Rob utelll993i However, for evo- 
lutionary paths starting with low initial obliquities, we can 
distinguish two different zones: one zone corresponding to 
slow initial rotation periods (Pi > 8 day) where the pro- 
grade rotation final state F$ is prevailing, and another zone 
for faster initial rotation periods (Pj < 8 day), where we 
find a mixture of the three attainable final states, Fq, F$ 
and F~ . To emphasize the chaotic behavior, we integrated 




Fig. 9. — Final states of Venus' spin with planetary perturbations 
for initial obliquity G [0°, 180°]) and period (Pi G [3d, 12 d]). 
For high initial obliquities, the final evolution of Venus remains 
essentially unchanged since none of the trajectories crossed the 
chaotic zone. The passage through the chaotic zone is reflected 
by the scattering of the final states in the left side of the picture. 
To emphasize the chaotic behavior, in the bottom left corner of 
picture (a), additional integrations were done with the same initial 
conditions, but with a differe nce of 10~ 9 in the initial ec centricity 
of Mars (b) and Neptune (c) dCorreia and Laskar 20033). 



twice more the zone with Pi < 8 day, with a difference of 
1CP 9 in the initial eccentricity of Mars (Fig.[9j3), and with 
a difference of 10~ 9 in the initial eccentricity of Neptune 
(Fig. [9};). The passage through the chaotic zone is reflected 
by the scattering of the final states in the left hand side of 
the picture. 

3.1.3. The Earth and Mars 

Contrary to Mercury and Venus, the Earth and Mars 
are not tidally evolved. For Mars, the tidal dissipation 
from the Sun is negligible. For the Earth, tidal dissi- 
pation is noticeable due to the presence of the Moon, 
but the Earth's spin is still far fr om the equilibrium (e.g. 
Neron de Surgy and Laskarlll997l) . Nevertheless, the spin 
axis of both planets is subjected to planetary perturbations 
and thus present some significant variations (Sect. 2.5). 

In the case of Mars, the presence of numerous secular 
resonances of the kind ijj = —Vf. (Eq.[79t induce large 
chaotic variations in the obliquity, which can evolve be- 
tween and 60 degrees (Fig.fTOb. At present, the obliquity 
of Mars is very similar to the obliquity of the Earth, which 
is a mere coincidence. Indeed, the obliquity of Mars has 
most certainly reached v alues larger than 45 degrees in the 
past ( Laskar et alJ2004a ). The high obliquity periods led to 
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20 40 60 80 100 
obliquity (deg) 

Fig. 10. — Frequency analysis of Mars' obliquity, (a) The fre- 
quency map is obtained by reporting in the ordinate the precession 
frequency value obtained for 1000 integrations over 56 Myr for the 
different values of the initial obliquity (abscissa). A large chaotic 
zone is visible, ranging from 0° to about 60°, with two distinct 
chaotic zones, Bl and B2. (b) Maximum and minimum values 
of the obliquity reached over 56 Myr. In (c), the power spectrum 
of the orbital forcing term (Eq.|76t is given in logarithmic scale, 
showing the correspondence of t he chaotic zone wi th the main sec- 
ular frequencies si, S2, S7, s$. dLaskar et alj2004al) . 



large climatic changes on Mars, with possible occurrence of 
large scale ice cycles where the polar caps are sublimated 
during high obliqui ty stages and the ice is deposited in the 
equat orial regions dLaskar et al.ll2002t iLevrard et al.1 12004 , 
2007bh . 

In the case of the Earth, the precession frequency is 
not in resonance with any orbital secular frequency (ip 
Vk). The obliquity of the Earth is then only subject to 
small oscillations of about 1.3 degrees around the mean 

periodicities around 40 000 
The small obliquity variations 
are nevertheless sufficiently important to induce substantial 
changes in the insolation received in summer in high lati- 
tude regions on the Earth, and the y are imprinted in the ge- 
ologi cal stratigraphic sequences dHavs et al.lll976t llmbrid 
1982b . 

Due to tidal dissipation in the Earth-Moon system, the 
Moon is moving a way from Earth at a 3.8 cm/yr rate 
(Dick ey et a 1. 1994), and the rotation rate of the Earth is 
slowing down (Eqs.[46] and [39l respectively). As a con- 
sequence, the torque exerted on the equatorial bulge of 
the Earth decreases and thus also the Earth precession fre- 
quency (Eqs.lTland fTSI) . Using the present dissipat i on pa- 
rameters of the Earth, iNeron de Surgy and Laskarl dl997t) 
found that after 1.5 Gyr, the spin of the Earth will enter 



value (23.3 degrees) with main ] 
years (|Laskar et alJl2004bh . The 
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Fig. 1 1 . — Example of possible evolution of the Earth's obliquity 
for 5 Gyr in the future, due to tidal dissipation in the Earth-Moon 
system. The background of the figure is obtained as the result of 
a stability analysis on about 250 000 numerical integrations of the 
obliquity of the Earth under planetary perturbations for 36 Myr, 
for various values of the initial obliquity of the Earth (x-axis) and 
various values of the precession constant (left y-axis). We ob- 
serve a very large chaotic zone (with stripes o n the figure) result- 
ing from overlap of orbita l secular resonances dLaskar etal J 19931 : 
lLaskar and Robutelll 19931) . The numerical integration is then con- 
ducted over 5 Gyr years in the future for the obliquity of the Earth, 
including tidal dissipative effects. The two bold curves correspond 
to the minimum and maximum values reached by the obliquity and 
the time scale is given in the right y-axis. As long as the orbits 
stays in the regular region, the motion suffers only small (and reg- 
ular) variations. As soon as the orbit enters the chaotic zone, very 
strong variations of the obliquity are observed, which wanders in 
all the c haotic zone, and very high values, close to 90 degrees, are 
reached dNeron de Surgy and Laskarll 19971) . 



a large chaotic zone of overlapping orbital secular reso- 
nances. From then, the Earth spin axis will evolve in a 
wildly chaotic way, with a possible range from to nearly 
90 degrees (Fig.QT,. 

The main difference between the Earth and Mars is thus 
due to the presence of the Moon, whose gravitational torque 
on the equatorial bulge of the Earth prevents the spin axis 
to evolve in a largely chaotic state. In absence of the Moon, 
the behavior of the spin axis of the Earth and Mars woul d 
be identical dLaskar et all 1 9931: lLaskar and Robute1ll993h . 

Depending on the orbital configuration of the exoplane- 
tary systems, we thus expect to find planets that would be 
either in a chaotic state as Mars or the moonless Earth, or 
in a regular state as the Earth with the Moon. It should 
be stressed, however, the presence of a large satellite is not 
mandatory in order to stabilize the spin axis. Since the sta- 
bility of the axis is very important for the exoplanet climate, 
planetary perturbations should be taken into consideration 
when searching for another Earth-like environments. 
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Table 3 

Constants for the Solar System outer planets. 



quantity Jupiter Saturn Uranus Neptune 



Po(h) 


9.92 


10.66 


17.24 


16.11 


P (g/cm 3 ) 


1.33 


0.69 


1.32 


1.64 


C/MR 2 


0.25 


0.21 


0.23 


0.24 


k 2 


0.49 


0.32 


0.36 


0.41 


Q (xlO 4 ) 


~3 


~2 


1~3 


1-30 



( Yoder 1995; Veeder et al. 1994; Dermott et al. 1988; Tittemore and Wisdom 1990; Banfield and Murray 1992) 



3.2. "Hot-Jupiters" 

One of the most surprising findings concerning exoplan- 
ets was the discovery of several giant planets with periods 
down to 3 days, tha t were designated by "Hot-Jupiters" (e.g. 



Santos et al.l 120051) . Many of the "Hot-Jupiters" were si- 



multaneously detected by transit method and radial Doppler 
shift, which allows the direct and accurate determination of 
both mass and radius of the exoplanet. Therefore, "Hot- 
Jupiters" are amongst the better characterized planets out- 
side the Solar System. 

Due to the proximity of the host star, "Hot-Jupiters" are 
almost certainly tidally evolved. Given the large mass of 
"Hot-Jupiters", they may essentially be composed of an ex- 
tensive Hydrogen atmosphere, similar to Jupiter or Saturn. 
As a consequence, despite the presence of an inner metallic 
core, "Hot-Jupiters" can be treated as fluid planets, and we 
may adopt a viscous model for the tidal dissipation (Eqs.l39l 
40b . Thermal atmospheric tides may also be p resent (e.g. 
Arras and Socratesl2 010, 2009 HGoodmanl2009l) . but we did 



not take thermal tides into account, as gravitational tides are 
so strong for "Hot-Jupiters", that they probably rule over all 
the remaining effects (see Sect. 3.3). 

3.2.1. Rotational synchronization 

The effect of gravitational tides over the obliquity is to 
straighten the spin axis (Eq.l45ll. so we will adopt e = 0° for 
the obliquity. Assuming that the eccentricity and the semi- 
major axis of the planet are constant, we can derive from 
expression i 



/2(e) ( ujq /2(e) 



Me) 



Me) 



exp(-t/T eq ) 



(93) 



where r^ 1 = Kf\(e)/G is the characteristic time-scale for 
fully despinning the planet. 

As for Mercury, the final equilibrium rotation driven 
by tides (t —> +00) is given by the equilibrium position 
Lu e /n = /2(e)//]. (e), which is different from synchronous 
rotation if the eccentricity is not zero (Fig.[4]i. Unlike Mer- 
cury, because "Hot-Jupiters" are assumed to be fluid, they 
should not present much irregularities in the internal struc- 
tures. Therefore (B — A)/C ~ and we do not expect 



"Hot-Jupiters" to be captured in spin-orbit resonances. In- 
deed, determination of second-degree harmonics of Jupiter 
and Satur n's gravity field from Pioneer a nd Voyager track- 
ing data ( Campbell and Anderson 1 19891) provided a crude 
estimate of the (B — A)/C value lower than — 10 -5 for 
Saturn and — 10~ 7 for Jupiter. The (B — A)/C values for 
Jupiter and Saturn are more than one order of magnitude 
smaller than the Moon's or Mercury's value, leading to in- 
significant chances of capture. Furthermore, the detection 
of an equatorial asymmetry is questionable. If an equatorial 
bulge originates from local mass inhomogeneities driven by 
convection, it is probably not permanent and must have a 
more negligible effect if averaged spatially and temporally. 

The time required for dampening the rotation of the 
planet depends on the dissipation factor k2 At (Eq.l43l. As- 
suming that "Hot-Jupiters" are similar to the Solar System 
giant planets, we can adopt k% = 0.4, and a range for Q 
from 10 4 to 10 5 (Table©. The Q-factor and the time lag 
At can be relied using expressions d22l i and ( 1341 : 



rAt . 



(94) 



Since we are using a viscous model, for which At is made 
constant, Q will be modified across the evolution as Q is 
inversely proportional to the tidal frequency a. The Q- 
factor for the Solar System gaseous planets is measured for 
their present rotation states, which correspond to less than 
one day (Table©. We may then assume that exoplanets 
should present identical Q values when they were rotating 
as rapidly as Jupiter, that is, Q^ 1 = uj^At. For Qo = 10 4 
and ojq — 27r/10h, we compute a constant At 0.57 s. 

We have plotted in Figure [12] all known exoplan- 
ets, taken from The Extrasolar Planets Encyclopedia 
( http://exoplanet.eu/), that could have been tidally evolved. 
We consider that exoplanets are fully evolved if their rota- 
tion rate, starting with an initial period of 10 h, is dampened 
toavalue such that |cj/n—/2(e)//i(e)| < 0.01. The curves 
represent the planets that are tidally evolved in a given time 
interval ranging from 0.001 Gyr to 10 Gyr . Figure [12] al- 
lows us to check whether the planet should be fully evolved. 
For a Solar type star, we can expect that all exoplanets that 
are above the 1 Gyr curve have already reached the equi- 
librium rotation ui e . On the other hand, exoplanets that are 
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semi-major axis (AU) 



Fig. 12. — Tidally evolved exospins with Qo = 10 4 and initial 
rotation period Po = 10 h. The labeled curves denote (in Gyr) 
the time needed by the rotation to reach the equilibrium (time- 
scales are linearly proportional to Qo)- We assumed Jupiter's 
geophysical parameters for all planets (Table[3} (updated from 
lLaskar and Corr eia 2004). 



below the 10 Gyr curve are probably not yet fully tidally 
evolved. As expected, all planets in circular orbits with 
a < 0.05 AU are tidally evolved. However, we are more 
interested on exoplanets further from the star with non zero 
eccentricity which are tidally evolved, since the rotation pe- 
riod is not synchronous, but given by expression d44i l. For 
instance, for the planet around HD 80606, the orbital pe- 
riod is 1 1 1.7 days, but since e = 0.92 we predict a rotation 
period of about 1.9 days. 

3.2.2. Cassini states 

Until now we have assumed th at the final obliquity of the 
planet is zero degrees. However. I Winn and Holmanl (12005) 
suggested that high obliquity values could be maintained 
if the planet has been trapped in a Cassini state resonance 
(Eq.|82t since the early despinning process. For small am- 
plitude variations of the eccentricity and inclination, the 
equilibrium positions for the obliquity are given by expres- 
sion ( f83l >. Unless a = \v).\, the state 1 is close to 0°, and 
state 3 is close to 180°. We thus focus only on state 2 that 
may maintain a significant obliquity. 

To test the possibility of capture in the high oblique 
Cassini state 2, we can consider a simple scenario where 
a "Hot-Jupiter" forms at a large orbital distance (~ sev- 
eral AUs) and migrates inward to the current position 
(~ 0.05 AU). Before the planet reaches typically ~ 0.5 AU, 
tidal effects do not affect the spin evolution, but the re- 
duction in the semi-major axis increases the precession 
constant (Eq.[T2l. so that the precession frequency ip may 
become resonant with some orbital frequencies (Eq.[80b. 
The passage throug h resonanc e generally causes the obliq- 
uity to change (Ward 1975; Ward and H amiltonl L2004; 
Hamilton and Wardl2004t iBoue et all 2009). rising the pos- 
sibility that the obliquity has a somewhat arbitrary value 



when the semi-major axis attains ~ 0.5 AU. 

Tidal effects become efficient for a < 0.5 AU and drive 
the obliquity to an equilibrium value cose w 2n(l+6e 2 )/w 
(Eq.l40T>. For initial fast rotation rates (u n), the equi- 
librium obliquity tends to 90°. As the rotation rate is de- 
creased by tides, the equilibrium obliquity is reduced to 
zero degrees (Eq.l45l). It is then possible that the obliquity 
crosses several resonances (one for each frequency Vk) in 
both ways (increasing and decreasing obliquity), and that 
a capture occurs. Inside the resonance island, the restor- 
ing torque causes the ob liquity to librate with amplitude 
dCorreia and Laskarll2003bh : 



COS £2 ± A COS £2 



1 



(95) 



Using a linear approximation of the tidal torque (Eq.l40l 
around the resonant obliquity £2, the probability of cap- 
ture in the Cassini state 2 can be estimated from the an- 
alytical approach fo r spin-orbit resonances (Eq.l73l, with 
dLevrard et al.ll2007al) 



Alo 
irVn 



(1 - 3 cos 2 £ 2 )^ + 2 cos £ 2 
7T sin 2 £2 (2 — cos £2 Lj/n) 



A cos £ 2 . (96) 



In FigureQ~3] we plotted the capture probabilities at 0.05, 
0. 1 and 0.5 AU as a function of the rotation period for dif- 
ferent amplitudes ( Jk) and frequencies (yk) ch aracteristic 
of the Solar System dLaskar and Robutellll993l) . As a rea- 
sonable example, we choose Vk = — 10"/yr, but the results 
are not affected by changes on this value. Assuming an ini- 
tial rotation period of 12 h, the capture is possible and even 
unavoidable if Jk > 0.1'Vyr at 0.5 AU. On the contrary, the 
chances of capture at 0.05 AU are negligible (< 1%) be- 
cause a decrease in the semi-major axis leads to an increase 
in the precession constant and reduces the width of the res- 
onance (Eq 



Theoretica l estimations can be compared 
with numerical simulations dLevrard et al. 2007a). To that 



purpose, the spin equations (Eqs. |79ll80t were integrated in 
the presence of tidal effects (Eqs. [391140b considering 1000 
initial precession angles equally distributed over — 2ir for 
each initial obliquity. Statistics of capture were found to be 
in good agreement with previous theoretical estimates. 

To test the influence of migration on the capture stability, 
additional numerical simulations were performed for vari- 
ous initial obliquities and secular perturbations over typi- 
cally ~ 5 x 10 7 yr. The migration process was simulated by 
an exponential decreasing of the semi-major axis towards 
0.05 AU with a 10 5 — 10 7 yr time scale. The obliquity libra- 
tions were found to be significantly shorter than spin-down 
and migration time scales so that the spin trajectory follows 
an "adiabatic invariant" in the phase space. Nevertheless, 
expression d40l ) indicates that the tidal torque dramatically 
increases both with spin-down and inward migration pro- 
cesses (de/dt oc a" 15 / 2 ^ 1 ). If the tidal torque exceeds 
the maximum possible restoring torque (Eq.|79l, the res- 
onant equilibrium is destroyed (the evolution is no longer 
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Fig. 13. — Obliquity capture probabilities in resonance under the effect of gravitational tides, as a functio n of the rotation perio d at a) 
0.5 AU b) 0.1 AU c) 0.05 AU. Jk is the amplitude of the secular orbital perturbations and Vk = — 10.0'Vyr dLevrard et alj2007ah . 



adiabatic). For a given semi-major axis, the stability con- 
dition requires then that the rotation rate must always be 
larger than a threshold value to cr it., which is always veri- 
fied if QJcrit. < /2 (e) /.fi (e). Th e stability condition can be 
simply written as (ILevrard et al ] l2007al) 



(97) 



tan(e) < J k x r eq , 

where r eq is the time scale of tidal despinning (Eq.1931. It 
then follows that the final obliquity of the planet cannot be 
too large, otherwise the planet would quit the resonance. 
For instance, taking J k = 1 "/yr at 0.05 AU (the highest 
value in Fig.Qj}, we need an obliquity e < 21°. For the 
more realistic amplitude J k =0.1 "/yr, the resonant obliq- 
uity drops to e < 2°. Such a low resonant obliquity at 
0.05 AU is highly unlikely, because, according to expres- 
sion (l83l l. the resonant state requires very high values of 
the orbital secular frequencies i\v k \ > 7.2 x 10 6 "/yr). At 
0.5 and 0.1 AU, critical obliquity values are respectively 
83° and 41° and require more reasonable orbital secular 
frequencies so that a stable capture is possible. In numer- 
ical simulation, the stability criteria for the final obliquity 
(Eq.l97b is em pirically retr ieved with an excellent agree- 
ment ( Levrard et alJl2007ah . When the obliquity leaves the 
resonance, the obliquity ultimately rapidly switches to the 
resonant stable Cassini state 1, which tends to 0° (Eq.[83l. 
We then conclude that locking a "Hot- Jupiter" in an oblique 
Cassini state seems to be a very unlikely scenario. 

3.2.3. Energy balance 

Tidal energy is dissipated in the planet at the expense 
of the rotational and orbital energy so that E = —CujCj — 
a(Gm+m) / '(2a 2 ). Replacing the equilibrium rotation given 
by equation (l44l in expression d46l > we obtain for the tidal 
energy: 



E = Kn l 



Me) 



fjje) 2x 2 
Me) 1 + 3 



(98) 
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Fig. 14. — Rate of tidal dissipation within HD 209458 b as a 
function of the eccentricity for 45° (solid thin line) and 90° (solid 
thick line) obliquity. The synchronous case is plotted with dashed 
lines for comparison . The dissipation factor Qo/fe is set to 10 6 
fcevrard et aU2007iJ) . 



or, at second order in eccentricity, 

Kn 2 



E 



1 + COS 2 £ 



[sin 2 e + e 2 



7 + 16 sin 2 e)] , (99) 



which is alwa ys larger than in the synchronous case (e.g. 
IWisdomll2004l) . In Figure[T4]the rate of tidal heating within 
a non-synchronous and synchronous planet as a function of 
the eccentricit y (0 < e < 0.25) is c ompared for two differ- 
ent obliquities (ILevrard et al . 2007a). The ratio between the 
tidal heating in the two situations is an increasing function 
of both eccentricity and obliquity. For e s» 0, as observed 
for "Hot-Jupiters", the ratio may reach ~ 1.3 and 2.0 at 
respectively 45° and 90° obliquity, not being significantly 
modified at larger eccentricity. 

We then conclude that planets in eccentric orbits and/or 
with high obliquity dissipate more energy than planets in 
synchronous circular orbits. This may explain why some 
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plane ts appear to be mo re inflated than initially expected 
(e.g. iKnutson et al]l2007l) . A correct tidal energy balance 
must then take into account the present spin and eccentricity 
of the orbit. 

3.2.4. Orbital circularization 

In Sect. 2.2.4 we saw that under tidal friction the spin of 
the planet attains an equilibrium position faster than the or- 
bit. As a consequence, we can use the expression of the 
equilibrium rotation rate (Eq.l44"li in the semi-major axis 
and eccentricity variations and find simplified expressions 
(Eqs.l5Tland[52"ll. Combining the two equations we get 



da 2e de 
~ = (1-e 2 ) 



whose solution is given by 



a f (l 



(100) 



(101) 



Replacing the above relation in expression (T52t we find a 
differential equation that rules the eccentricity evolution: 



e = -K f 6 (e)(l-e 2 fe, 
where Ko is a constant parameter: 

21k 2 GmlR 5 



K = At- 



The solution of the above equation is given by 

F(e) = F{e )e- K ^ , 



(102) 



(103) 



(104) 



where F(e) is an implicit function of e, which converges to 
zero as t — > +oo. For small eccentricities, we can neglect 
terms in e 4 and F(e) = e\7 — 9e 2 1 ~ x / 2 . The characteristic 
time-scale for fully dampening the eccentricity of the orbit 
is then T or b ~ 1/Kq, and the ratio between the spin and 
orbital time-scales 



' eq 
7~orb 



(105) 



Since aj — a(l — e 2 ), for initial very eccentric orbits the 
two time-scales become comparable. 

We have plotted in Figure Q3] all known exoplan- 
ets, taken from The Extrasolar Planets Encyclopedia 
(http://exoplanet.eu/), whose orbits could have been 
tidally evolved. We consider that they are fully evolved 
if the eccentricity is dampened to a value e < 0.01. The 
curves represent the time needed to damp the eccentric- 
ity starting with the present orbital parameters, for time 
intervals ranging from 0.001 Gyr to 100 Gyr. Figure [T31 
allows us to simultaneously check whether the planet is 
tidally evolved, and the time needed to fully damp the 
present eccentricity. All planets in eccentric orbits expe- 
rience stronger tidal effects because the planet is close to 
the star at the periapse. As a consequence, tidal friction 
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Fig. 15. — Tidally evolved orbits of exoplanets with Qo = 10 4 . 
The labeled curves denote (in Gyr) the time needed to circularize 
the orbits of the planets (e < 0.01) (time-scales are linearly pro- 
portional to Qo). We assumed Jupiter's geophysical parameters 
for all planets (Tabled. 



can, under certain conditions, be an important mechanism 
for the formation of "Hot-Jupiters" (see Sect. 3.2.5). 

According to Figure[l5] a significant fraction of exoplan- 
ets that are close to the host star (a < 0.1 AU) still present 
eccentricities up to 0.4, although tidal effects should have 
already damped the eccentricity to zero. Observational er- 
rors and/or weaker tidal dissipation (Qo 3> 10 4 ) can be 
a possible explanation, but they will hardly justify all the 
observed situations. A more plausible explanation is that 
the eccentricity of the exoplanet is being excited by grav- 
itational perturbations from an outer planetary companion 
(Sect. 2.5). Indeed, the eccentricity of a inner short-period 
planet can be excited as long as its (non-resonant) outer 



Mardlingl(l2007h has 



companion s eccentricity is non-zero, 
shown that the eccentricity of the outer planet will decay 
on a time-scale which depends on the structure of the in- 
ner planet, and that the eccentricities of both planets are 
damped at the same rate, controlled by the outer planet 
(Fig.fToT). The mechanism is so efficient that the outer planet 
may be an Earth-mass planet in the ''habitable zones" of 
some stars. As a consequence, the evolution time-scale for 
both eccentricities can attain the Gyr instead of Myr, which 
could explain the current observations of non-zero eccen- 
tricity for some "Hot-Jupiters". 

3.2.5. Kozai migration 

In current theories of planetary formation, the region 
within 0.1 AU of a protostar is too hot and rarefied for a 
Jupiter-mass planet to form, so "Hot-Jupiters" likely form 
further away and then migrate inward. A significant fraction 
of " Hot-Jupiters" has been fou nd in systems of binary stars 
(e.g. lEggenberger et al.1 120041) . suggesting that the stellar 
companion may play an important role in the shrinkage of 
the planetary orbits. In addition, close binary star systems 
(separation comparable to the stellar radius) are also often 
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years years 

Fig. 16. — Tidal evolution of the eccentricities of planet HD 209458 b (e p ) perturbed by a 0.1Mj up companion at 0.4 AU (e c ). The 
dissipation factor Q v = 100 of the observed planet is set artificially low in order to illustrate the damping process (time-scales are 
linearly proportional to Q p ). Top figures show the first stages of the evolution. The change in grayscale shows the transition of the 
eccentricity from the circulation phase to the libration phase. The diamond represents the moment where the eccentricity librations are 
damped. Bottom figures show the final evolution of the eccentricities. The p resence of a co mpanion result that e v decays at the same 
rate of e c , while the dissipation rate for e c is controlled by Q v and not by Q c |Mardling 20o3). 



accom panied by a third star. For instance, iTokovinin et al 
(2006) bund that 96% of a sample of spectroscopic bina- 
ries with periods less than 3 days has a tertiary component. 
Indeed, in some circumstances the distant companion en- 
hances tidal interactions in the inner binary, causing the bi- 
nary orbital period to shrink to the currently observed val- 
ues. 

Three-body systems can be stable for long-time scales 
provided that the system is hierarchical, that is, if the system 
is formed by an inner binary (star and planet) in a nearly 
Keplerian orbit with a semi-major axis a, and a outer star 
also in a nearly Keplerian orbit about the center of mass of 
the inner system with semi-major a' 3> a. An additional 
requirement is that the eccentricity e' of the outer orbit is 
not too large, in order to prevent close encounters with the 
inner system. In this situation, perturbations on the inner 
planetary orbit are weak, but can have important long-term 



effects (Chapter 10: Non-Keplerian Dynamics). 

The mos t striking effect is know n as the Lidov-Kozai 
mechanism ( Kozaill962 ; Lidovl 19621) . which allows the in- 
ner orbit to periodically exchange eccentricity with incli- 
nation. Even at large distances (a! >1000 AU), the outer 
star can significantly perturb the planetary orbit as long as 
the two orbital planes are initially inclined to each other 
more than I > 39.2° (that is, for cos / < (3/5) 1 / 2 ). When 
/ < 39.2° there is little variation in the planet's inclination 
and eccentricity. Secular effects of the Lidov-Kozai type 
can then produce large cyclic variations in the planet's ec- 
centricity e as a result of angular momentum exchange with 
the companion orbit. Since the z-component of the planet's 
angular momentum must be conserved and a is not modi- 
fied by the secular perturbations, the Kozai integral 



L K = (l-e 2 ) 1/2 cos/ 



(106) 
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Fig. 17. — Possible evolution of the planet HD 80606b initially in an orbit with a = 5 AU, e = 0.1, and I = 85.6°. The stellar 
companion is supposed to be a Sun-like star at a' = 1000 AU, and e' = 0.5. Th e diamonds mark the current position of HD 80606 b 
along this possible evolution dWu and Murr ay 2003; Fabrvckv and Tremainell2007h . 



is co nserved during the oscillations ( Lidov and Ziglinl 
1976h . Maxima in e occur with minima in /, and vice 
versa. If the inner orbit is initially circular, the maximum 
eccentricity achieved in a Kozai cycle is e max = (1 — 
(5/3) cos 2 I) 1 / 2 and the oscillation period of a cycle is ap- 
proximately P' 2 /P (Kiseleva et al. 1998). The maximum 
eccentricity of the inner orbit in the Kozai cycle will remain 
fixed for different masses and distances of the outer star, but 
the period of the Kozai cycle will grow with a' 3 . Kozai cy- 
cles persist as long as the perturbation from the outer star is 
the dominant cause of periapse precession in the planetary 
orbit. However, small additional sources of periapse pre- 
cession such as the quadrupole moments, additional com- 
panions, general relativity or even tides can compensate the 
Kozai precessio n and suppress the eccentricity/inclin ation 
oscillations (e.g. Migaszewski and Gozdziewski 2009). 

Because the Lidov-Kozai mechanism is able to induce 
large eccentricity excitations, a planet in an initial almost 
circular orbit (for instance a Jupiter-like planet at 5 AU 
around a Sun-like star) can experiment close approaches to 
the host star at the periapse when the eccentricity increases 



to very high values. As a consequence, tidal effects increase 
several orders of magnitude and according to expression 
(l5lT i the semi-major axis of the orbit will decrease and the 
planet migrate inward. At some point of the evolution, the 
periapse precession will be dominated by other effects and 
the eccentricity oscillations suppressed. From that moment 
on, the eccentricity is damped according to expression d52| > 
and the final semi-ma jor axis given by af = a(l — e 2 ) 



Ford and Rasiol d2006l) have derived that tidal evolution of 



high eccentric orbits would end at a semi major axis a/ 
equal to about twice the Roche limit Rl- Indeed, at the 
closest periapse distance, attained for e s» 1, we will have 
a(l — e) = Rl, and thus at = (1 + b)Rl ~ 2i?£. 

In Figure[l7]we plot an example of combined Kozai-tidal 
migration of the planet HD 80606 b. The planet is initially 
set in an orbit with a = 5 AU, e = 0.1, and I = 85.6°. 
The stellar companion is supposed to be a Sun-like star 
at a' = 1000 AU, and e' = 0.5 dWu and Murravl 12003 



Fabrvckv and T remaine 2007). Prominent eccentricity os 



cillations are seen from the very beginning and the energy 
in the planet's spin is transferred to the orbit increasing the 
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semi-major axis for the first 0.1 Gyr (Eg. [391). As the equi- 
librium rotation is achieved (Eq.l44l> the orbital evolution 
is essentially controlled by equations (l51~l l and (l52l . whose 
contributions are enhanced when the eccentricity reaches 
high values. The semi-major axis evolution is executed by 
apparent "discontinuous" transitions precisely because the 
tidal dissipation is only efficient during periods of high ec- 
centricity. As dissipation shrinks the semi-major axis, peri- 
apse precession becomes gradually dominated by relativity 
rather than by the third body, and the periapse starts cir- 
culating as the eccentricity passes close to at 0.7 Gyr. 
Tidal evolution stops when the orbit is completely circular- 
ized. The present semi-major axis and eccentricity of planet 
HD 80606b are a = 0.45 AU and e = 0.92, respectively, 
meaning that the tidal evolution on HD 80606 b is still un- 
der way (Fig.lT3Tl. The final semi-major axis is estimated to 
about a/ = 0.07 AU, which corresponds to a regular "Hot- 
Jupiter". 

3.3. "Super-Earths" 

After a significant number of discoveries of gaseous 
giant exoplanets, a new barrier has been passed with the 
detections of several exoplanets in the Neptune and even 
Earth-mass (M m ) regime: 2 - 12 jRivera et dJl2005t 
Lovis et all 120061 lUdry etalJ l2007t iBonfils et all 120071 . 



that are commonly designated by "Super-Earths". If the 
commonly accepted core-accretion model can account for 
the formation of "Super-Earths", resulting in a mainly 
icy/rocky composition, the fraction of the residual He-H2 
atmospheric envelope accreted during the planet migration 
is not tight ly constrained for p lanets more massive than the 
Earth (e.g. lAlibert et al . 2006). A minimum mass of below 
10 is usually considered t o be the bounda ry between 
terrestrial and giant planets, but|RafikoyJ (120061) found that 
planets more massive than 6 M ffi could have retained more 
than 1 of the He-H2 gaseous envelope. For compari- 
son, masses of Earth's and Venus' atmosphere are respec- 
tively ~ 10~ 6 and 10~ 4 times the planet's mass. Despite 
significant uncertainties, the discoveries of "Super-Earths" 
provide an opportunity to test some properties that could be 
similar to those of the more familiar terrestrial planets of 
the Solar System. 

Because some of t he "Super-Earths" are potentially i n 
the "habitable zone" dUdrv et alj|2007l: ISelsis et alj|2007l) . 
the present spin state is an important factor to constrain the 
climates. As for Venus, thermal atmospheric tides may have 
a profound influence on the spin of "Super-Earths". How- 
ever, the small eccentricity approximation done for Venus 
(Eq.l87t may no longer be adequate for "Super-Earths", 
which exhibit a wide range of eccentricities, orbital dis- 
tances, or central star types. Although our knowledge of 
"Super-Earths" is restricted to their orbital parameters and 
minimum masses, we can attempt to place new constraints 
on the surface rotation rate, assuming that "Super-Earths" 
have a dense atmosphere. 

As for Venus, the combined effect of tides is to set the 





a>2 








-v.. a>2 


~cos/n= 


1.92 




r 







2 -1 





1 2 


3 


















a)s/n=0.4( 


) 






0.5 


1 1.5 





6 
C 


= 0.0 

= 0.1 
= 0.2 


a)s/n=0.0f 


i 





1 

co/n 



Fig. 18. — Evolution of Cj (Eq. ll07t with uj/n for different at- 
mospheric strengths (uj s /n = 1.92, 0.40, 0.05) and eccentricities 
(e = 0.0, 0.1, 0.2). The top picture with e = is the same as Fig- 
ure[8]for Venus. The equilibrium rotation rates are given by to = 
and the arrows indicate whether the equilibrium position is stable 
or unstable. For u s /n > 1, we have two equilibrium possibili- 
ties, wj, one of which corresponds to a retrograde rotation (as for 
Venus). For u) a /n < 1, retrograde states are not possible, but we 
can still observe final rotation rates < n. For eccentric orbits, 
because of the harmonics in a = 2ui — n and a = 2oj — 3n, we 
may have at most four different final possibilities (Eq. ll 10t . When 
uj s /n becomes extremely small, which is the case for the present 
observed exoplanets with some ec centricity (Tab led, a single fi- 
nal equilibrium is possible for uj^ dCorreia et al . 2008). 



final obliquity at 0° or 180° dCorreia et alj|2003l) . Adopting 
a viscous dissipation model for tidal effects ( Eq.[37b a nd the 
"heat ing at the ground' model dDobrovolskis and Ingersoll 
1980) for surface pressure variations (Eq.|68l), the average 
evolution of the rotation rate is then obtained by adding the 
effects of both tidal torques acting on the planet. From ex- 
pressions d23l ) and d60b we get for e — 0° and to the second 
order in the eccentricity: 

— — = ui — (l + 6e 2 ) n — uj s [(l - ^-e 2 ) sign(w - n) 

- e 2 sign(2w - n) + 9 e 2 sign(2w - 3n)] , (107) 



where 



F s K a At a L* R 

oc a . 

m+ m 



\&H k 2 K g M g " il0S ' 

When e = 0, we saw in the case of Venus that final 
positions of the rotation rate at zero obliquity are given by 
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Fig. 19. — Equilibrium positions of the rotation rate for "Super-Earths" as a function of the product am* for three different values of 
the eccentricity (e = 0.0, 0.1, 0.2). Each curve corresponds to a different final state (dotted lines for cuf and solid lines for (jJ^ 1 ). For 
e ~ (case of Venus), we always count two final states that are symmetrical about n. For small values of am*, the two equilibrium 
possibilities are so close to n that the most likely scenario for the planet is to be captured in the synchronous resonance (case of GJ 58 1 e). 
As we increase the eccentricity, we can count at most three final equilibrium rotations, depending on the value of ui s /n (computed from 
Eq. ll 14t . When e ~ 0.2, only one equilibrium state exists for am* < 0.3, resulting from to B /n < 6e 2 (l + e 2 /2). This is the present 
situation of \i Araec and most of the "Super-Earths" listed in Tablel4l(Corr eia et alj2008h . 



(Eq®: 



ui 



(109) 



i.e., there are two final possibilities for the equilibrium ro- 
tation of the planet, given by cj ± = n ± uj s . When e ^ 0, 
the above expression (U09t is no longer valid and additional 
equilibrium positions for the rotation rate may occur. For 
moderate values of the eccentricity, from expression ( 1107t 
we have that the effect of the eccentricity is eventually to 
split each previous equilibrium rotation rate into two new 
equilibrium values. Thus, four final equilibrium positions 
for the rotation rate are possible (eight if we consider the 
case e = 180°), obtained with Co = (Fig.[T8Ti: 



w i 2 = n± ui s 



'5 ± 
°1 ,: 



(HO) 



The Earth and Venus are the only planets that can be 
included in the category of "Super-Earths" for which the 
atmosphere and spin are known. Only Venus is tidally 
evolved and therefore suitable for applying the above ex- 
pressions for tidal equilibrium. We can nevertheless in- 
vestigate the final equilibrium rotation states of the already 
detected "Super-Earths". For that purpose, we considered 
only exoplanets with masses smaller than 12 Af® that we 
classified as rocky planets with a dense atmosphere, al- 
though we stress that this mass boundary is quite arbitrarily. 

Using the empirical mass-luminosity relation L* tx Aff 
(e.g. lCester et al.ll983l) and the mass-radius relatio nship for 



terrestrial planets R oc m 274 dSotin et al.ll2007l) . expres- 
sion d 1 08t can be written as: 



with 



and 



1 , 41 
S 1 = 6n+-u; s , = 6n — co s , (111) 



5 5 
62=6n+-uj s , 62=6n--u s . (112) 

Because the set of Lof 2 values must verify the additional 
condition 



< n/2 < U) x < n < lu x < in/2 < u\ 



(113) 



the four equilibrium rotation states cannot, in general, ex- 
ist simultaneously, depending on the values of co s and e. In 
particular, the final states lo^ and io^ can never coexist with 
lo 2 . At most three different equilibrium states are therefore 
possible, obtained when io s /n is close to 1/2, or more pre- 
cisely, when 1/2 - 19e 2 /4 < u s /n < 1/2 + 17e 2 /4. 
Conversely, we find that one single final state ui x = (1 + 
6e 2 ) n+(l-41e 2 /2)w s exists when ui s /n < 6e 2 (l+e 2 /2). 



uj s /n = k (am*) 



2.5 m -0.726 



(114) 



where k is a proportionality coefficient that contains all the 
constant parameters, but also the parameters that we are un- 
able to constrain such as Hq, ki, At g or At a . In this con- 
text, as a first order approximation we consider that for all 
"Super-Earths" the parameter n has the same value as for 
Venus. Assuming that the rotation of Venus is presently sta- 
bilized in the oj~ final state, that is, 2-k — —243 days 
(Car penter 1970i) . we compute 2k/uj s — 116.7days. Re- 
placing the present rotation in expression (II 14b . we find for 



Venus that k = 3.723 M, 



0.720 



Mq 2 - 5 AU 



-2.5 



We can then 



estimate the ratio ui s /n for all considered "Super-Earths" 
in order to derive their respective equilibrium rotation rates 
(Table!). 

The number and values of the allowed equilibrium ro- 
tation states are plotted as a function of aM* for different 
eccentricities in Figure[l9] All eccentric planets have a ratio 
uj s /n that is lower than 6 x 10~ 3 (Table|4]i, which verifies 
the condition oj s /n < 6e 2 (l + e 2 /2). As a consequence, 
only one single final state exists uif jn k, (1 + 6e 2 ), cor- 
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Table 4 

Characteristics and equilibrium rotations of some "Super-Earths" with masses lower than 12 M ( 



Name m* Age r eq msini a 

[M Q ] [Gyr] [Gyr] [Af ffi ] [AU] 



2n/n 
[day] 



2ty/ui 2 
[day] 



[day] 



2?r/w+ 
[day] 



[day] 



Venus 


1 nn 
l.UU 


1 c 
4. J 


Z.J 


U.oz 


n t)q 
KJ./Zj 


U.UVI 


L.yZ 


Zz4. / 


iai n 








1 nn 
l.UU 


4. J 


1 O 


1 nn 
1 .UU 


1 nnn 
l.UUU 


U.U1 / 


D. 1 J 


juj.j 


1 11 Q 






^JSSle 


0.31 


7-11 


io- 7 


1.94 


0.03 





10" 5 


3.4087 




3.4088 


3.4087 


2 HD 40307 b 


0.77 




io- 7 


4.2 


0.047 





0.0003 


4.2413 




4.2427 


4.240 


^JSSlc 


0.31 


7-11 


10~ 5 


5.36 


0.07 


0.17 


io- 4 


12.14 






10.6335 


3 GJ876d 


0.32 


9.9 


IO -8 


6.3 


0.021 


0.14 


10" 6 


1.9649 






1.7822 


2 HD 40307 c 


0.77 




10~ 5 


6.9 


0.081 





0.0009 


9.5956 




9.6042 


9.5871 


^JSSld 


0.31 


7-11 


0.02 


7.09 


0.22 


0.38 


0.0011 


67.6918 






47.8226 


4 HD 181433 b 


0.78 




io- 5 


7.5 


0.08 


0.396 


0.0008 


9.3579 






6.535 


5 GJ176b 


0.5 




10~ 5 


8.4 


0.066 





0.0002 


8.7583 




8.7596 


8.7568 


2 HD 40307 d 


0.77 




io~ 4 


9.2 


0.134 





0.0025 


20.4175 




20.4696 


20.3656 


6 HD 7924 b 


0.83 




10~ 6 


9.26 


0.057 


0.17 


0.0004 


5.4493 






4.7688 


7 HD 69830 b 


0.86 


4-10 


io- 5 


10.2 


0.079 


0.10 


0.0008 


8.6625 






8.1995 


8 /i Araec 


1.1 


6.41 


io- 5 


10.6 


0.091 


0.172 


0.0021 


9.5505 






8.13313 


9 55Cnce 


1.03 


5.5 


io- 7 


10.8 


0.038 


0.07 


0.0002 


2.6659 






2.592 


10 GJ674b 


0.35 


0.1-1 


IO" 7 


11.09 


0.039 


0.2 


10~ 5 


4.7549 






4.0138 


7 HD 69830 c 


0.86 


4-10 


10~ 3 


11.8 


0.186 


0.13 


0.0064 


31.5943 






28.8691 



76.8 
77.1 



(*) Moon tidal effec ts w ere not included. r eo was computed wit h k 



iMavor et alj d2009ah: r2llMavor et al.1 d2009bh: mlCorreia et al 



Howard et aljj2009h ; r7l lLovis et alj<2006b : r8l |Pepe et al. 



(2010); [4] 



(2007); [9] 



1/3 and At a = 640 s (E arth 's values). Refe r ences : [1] 
iBouchv et al.U2009h ; r5l lForveille etai](l2009h ; [6] 



Fischer et alH2008h ; [10] Bonfil s et alj I ( I2007h 



responding to the equilibrium rotation resulting from gravi- 
tational tides (Eq.l44"l>. The main reason is that the effect of 
atmospheric tides is clearly disfavored relative to the effect 
of gravitational tides on "Super-Earths" discovered orbiting 
M-dwarf stars: the short orbital periods strengthens the ef- 
fect of gravitational tides, which are proportional to 1 /a 6 , 
while the effect of thermal tides varies as 1/a 5 . Moreover, 
the small mass of the central star also strongly affects the 
luminosity received by the planet and hence the size of the 
atmospheric bulge driven by thermal contrasts. 

For the planets with nearly zero eccentricity (GJ581 e, 
HD 40307 b,c,d, and GJ176b), two equilibrium rotation 
states ujf are possible. However, the two final states cof 1 
are so close to the mean motion n, that the quadrupole mo- 
ment of inertia (B — A)/C will probably capture the rota- 
tion of the planet in the synchronous resonance. We then 
conclude that "Super-Earths" orbiting close to their host 
stars (in particular M-dwarfs), will be dominated by grav- 
itational tides and present a final equilibrium rotation rate 
given by oj e /n s» f% (e) //i (e) (Fig.|Ul, or present spin-orbit 
resonances like Mercury. 

4. FUTURE PROSPECTS 

Th e classical theory of tides initiated by Darwinl (1880, 
1 19081) is sufficient to understand the main effects of tidal 
friction upon planetary evolution. However, the exact 
mechanism on how tidal energy is dissipated within the in- 



ternal laye rs of the planet remains a challenge for planetary 



scientists. iKaulal (11964) derived a generalization of Dar- 
win's work, with consideration of higher order tides and 
without the adoption of any dissipation model. The tidal 
potential is described using infinite series in eccentricity 
and inclination, which is not practical and can only be cor- 
rectly handled by computers. Ever since many efforts have 
been done in order to either simplify the tidal equations, 
or to correctly model the tidal dissipation (for a review see 



Ferraz-Mello et al . 2008; Efroimsky and Wimamsll2009l) . 

Many Solar System phenomena have been successfully 
explained using the existent tidal models, so we expect that 
they are suitable to describe the tidal evolution of exoplan- 
ets. Nevertheless, many exoplanets are totally different 
from the Solar System cases, and we cannot exclude to ob- 
serve some unexpected behaviors. For instance, it is likely 
that dissipation withi n "Hot-Jupi ters" is closer to dissipa- 
tion within stars (e.g. lZahnll 1975b . while dissipation within 
"Super-Eart hs" is closer to dissip ation observed for rocky 



planets (e.g. iHenning et al. 2009). It is then necessary to 



continue improving tidal models in order to get a more re- 
alistic description for each planetary system. In particular, 
a correct description of the tidal dissipation and on how it 
evolves with the tidal frequency is critical for the evolution 
time-scale. 

The orbital architecture of exoplanetary systems is rela- 
tively well determined from the present observational tech- 
niques. However, the spins of the exoplanets are not easy to 
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measure, as the light curve coming from the planet is always 
dimmed by the star light. The continuous improvements 
that have been made in photometry and spectrography let 
us believe that the determination of exoplanets' spins can 
be a true possibility in the near future. In particular, infra- 
red spectrographs are being developed, which will allow to 
acquire spectra of the pl anets if we manag e to subtract the 
stellar contribution (e.g. lBarnes et al 

Some additional methods for detecting the rotation 
and/or the obliquity of exoplanets have also been tested 
and suggested so far. For instance, indirect sensing of 
the planetary gravitational quadrupole and shap e, which is 
linked to both spin rate and obliq uity ("e.g. ISeager and Huil 
20021 iRagozzine and Wold 12009). or transient heating of 
one face of the planet, which then spins into and out of 
view, as it has been attempted for the system HD 80606 
(Laughlin et al. 2009). The effect of planetary rotation on 
the transit spectrum of a giant exoplanet is another possibil- 
ity. During ingress and egress, absorption features arising 
from the planet's atmosphere are Doppler shifted by of or- 
der the planet's rotational velocity (^l-2kms -1 ) relative 
to where they would be if the planet were not rotating (e.g. 



Spieg el et alJ l2007). Finally, for planets whose light is spa- 
tially separated from the star, variations may be discernible 
in the light curve obtained by low -precision photometry due 
to meteor ological variabi lity, composition of the surface, or 
spots (e.g. lFordetai1l200ll) . 

Although the spin states of exoplanets cannot be mea- 
sured, for exoplanets that are tidally evolved we can still 
try to make predictions for the rotation rates. When the ec- 
centricity is large, the rotation of many of the observed exo- 
planets can still be tidally evolved even if the planets are not 
very close to their central stars (Fig. [12]). For tidally evolved 
"Hot-Jupiters", we can conjecture that the rotation periods 
are the limit values P or i, x /i(e)//2(e) (Fig.|4j». It becomes 
a new challenge for the observers to be able to confirm these 
predictions. 

Thermal atmospheric tides may very well destabilize the 
tidal equilibrium from gravitational tides and create addi- 
tional possible stable limit values, with the possibility of 
retrograde rotations, as for planet Venus (FigH). Thermal 
tides should be particularly important for "Super-Earths", 
which are expected to have a distinct rocky body sur- 
rounded by a dense atmosphere. In a paradoxical way, the 
final rotation rate of "Super-Earths" are the most difficult 
to predict, as the equilibrium configurations depend on the 
composition of the atmospheres. Thermal tides are never- 
theless more relevant for exoplanets that orbit Sun-like stars 
at not very close distances, like Venus (FigfT9ll. 

We also assumed that the final obliquity of exoplanets is 
either 0° or 180°, as the two values represent the final out- 
come of tidal evolution. However, each planetary system 
has its own architecture, and planetary perturbations on the 
spin can lead to resonant capture in a high oblique Cassini 
states or even to chaotic motion. Thus, the final spin evolu- 
tion of a planet cannot be dissociated from its environment, 
and a more realistic description of exoplanets rotation can 



only be achieved with the full knowledge of the system or- 
bital dynamics. 
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